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Abstract 

We study the phase structure and Hall conductance quantization in weakly 
coupled multi-layer electron systems in the integer quantum Hall regime. We 
derive an effective field theory and perform a two-loop renormalization group 
calculation. It is shown that (i) finite interlayer tunnelings (however small) 
give rise to successive metallic and insulating phases and metal-insulator tran- 
sitions in the unitary universality class, (ii) The Hall conductivity is not 
renormalized in the metallic phases in the 3D regime, (iii) The Hall con- 
ductances are quantized in the insulating phases. In the bulk quantum Hall 

phases, the effective field theory describes the transport on the surface. 
PACS numbers: 73.50.Jt, 73.40.Hm, 72.15.Rn, 71.30.+h 



1 



Typeset using REVR3X 



The quantum Hall effect (QHE) in a two-dimensional electron gas (2DEG) has led to 
many new physical concepts and principles M. The main part of the phenomenology is 
that in strong magnetic fields, a 2DEG exhibits continuous zero-temperature phase transi- 
tions between successive quantum Hall states of vanishing dissipation and quantized Hall 
resistances. It is interesting to ask what happens to the physics associated with the QHE 
in dimensions greater than two 0. Experimentally, two classes of quasi-three-dimensional 
materials have been found to show integer-quantized Hall plateaus: multi-layer quantum 
wells or superlattices formed by GaAs/AlGaAs graded heterostructures fUJiJ] and molecular 
crystals (TMTSF^X (X=PF6,C104) ||. In this paper, we concentrate on the former which 
is a natural generalization of the QHE above two dimensions. Thus, the changes in the 
phase structure and the properties of the phase transitions in quantum Hall layers coupled 
by weak interlayer tunnelings are the concerns of the present paper. Moreover, we focus on 
the integer quantum Hall regime and ignore the effects of electron-electron interactions Q. 

We shall follow the approach of Chalker and Dohmen who generalized the Chalker- 
Coddington network model ]7J for the integer QHE (IQHE) in a 2DEG to layers of networks 
coupled by interlayer tunnelings ||. The advantage of this approach is that the single 
layer network model is known to correctly describe the universality class of the 2D integer 
quantum Hall transitions. Chalker and Dohmen performed numerical transfer matrix cal- 
culations and demonstrated the existence of three phases: insulator, metal, and quantized 
Hall conductor, and extended surface states in the quantized Hall state. In this paper we 
provide an analytical treatment using the effective field theory representation of the network 
model IP-pf. We first show that the long wavelength transport properties are governed by 
a 3D anisotropic unitary nonlinear a-model (NLcrM). The anisotropic couplings are the dis- 
sipative conductivities, whereas the Hall conductivity enters as a coupling to the layered 
sum of the 2D topological term. The renormalization group (RG) is then used to study the 
crossover between two and three dimensions. In the 3D regime, the RG flow equations for 
the conductances are calculated to two-loop order to determine the phase structure in the 
plane spanned by the Fermi energy and the interlayer tunneling. The results show that a 
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finite interlayer coupling (however weak) leads to metallic and insulating phases and metal- 
insulator transitions in the unitary universality class. Furthermore, the Hall conductivity 
is found to be unrenormalized by localization effects in the 3D regime. We show that the 
Hall conductance is quantized in the insulating phases provided that the above results hold 
to all orders in the RG. Finally, we demonstrate that in the quantum Hall phases, the field 
theory reduces to the one appropriate for the coupled edge states on the surface. 

Following Ref. ||10|| , the Hamiltonian for N layered networks in the (x,y)-plane coupled 
in the z-direction by interlayer tunneling t± is, 

H = XX -1 ) / d y^(x, y, z) [id y - V x , y>z ] ip(x, y, z) 

x,z J 

-53 / d v ty( x + h y, z)^(x, y, z ) + h - c - 

x,z 

-5Z / d y t± [i>\x, y,z + 1)^0, y, z) + h.c.] . (1) 

x,z J 

Here ip^ creates an electron traversing the edges of the Hall droplets as modeled by the 
links of the network. t x = t[l — 8(— 1) E ], where 5 measures the distance of the Fermi energy 
(Ep) relative to the center of the Landau level {E c ), represents the quantum tunneling 
amplitudes at the saddle points of the random potential, i.e. at the nodes of the network. V 



is a local random variable that generates the link Aharonov-Bohm phases, < V x ^^ z V x i t 



> 



U8 XiX >8 ZjZ >8(y - y'). 

For t± = 0, Eq. ([[]) describes N decoupled 2D networks, each undergoes a quantum 
Hall transition as 5 is varied. In this case, quench averaging over V and regarding y as the 
Euclidean time r, it has been shown that the original network model corresponds to a half- 
filled ID U(2n) Hubbard model in the limit n — > ||10|| . The 2D quantum Hall transition 
is then equivalent to the dimerization transition of the Hubbard chain EDI. Generalizing to 
t± 7^ 0, we obtain the generating functional Z = JV[^>, ip) exp[J dr X^z^'Sp — Hq^, ip))] 
in the form of a 2+1-dimensional Euclidean action if we let ip p — > ip p (iip p ) and ip p — > —i^ p {^ p ) 
for even (odd) x [I(J, 



S 



dr 



2 ^ad T ^a + H\ 



(2) 



Here 77 is a positive infinitesimal, a = (a,p) are the replica index a = 1, . . . ,n and energy 
index p = +(— ) for the advanced (retarded) channels, and S p = sgn(p). The resulting 
Hamiltonian H in Eq. (|2[) corresponds to an interacting quantum theory in two spatial 
dimensions, 

U 2 

x,z x,z 
+ £ tt^-lf tyl(x, Z + l)^*, *) + h.C.] - T) Y,(-l) X S p ^a, (3) 

£,2 X,Z 

where sums over repeated indices are implied. Note that Eq. (|3|) is not the usual quasi- ID 
U(2n) Hubbard model for the form of the interchain couplings. 

We now derive the effective low energy, long wavelength theory. The partition function 
for Eq. ([5]) can be written as Z = JV[ip , ip]V[A] exp(— S), 

s= f dr^2-v(-^ x s P ^ a + s + s I + s L 

J x,z 

S = dT Yl [i>ad T lpa ~ t x {i> a {x + l)lp a (x) + C.C. 
J x,z 

Si= jdrY, 

J x,z 

S± = [dr^2it L (-l) x Mz + l)Mz) + cc 



2U 



i) a A ab 4> b - -A ab 5 ab 



(4) 



where A(x, r, z) is a matrix Hubbard-Stratonovich field. As usual, at a mean-field level 
A°& = U(4> b 4> a - 6^/2), which is easily solved to give A° ab (x) = A (-l) x A ab with A ab = 
S p 5 ab . The massless fluctuations beyond the mean- field theory can be represented by slowly- 
varying unitary rotations of the "staggered magnetization" A® b . Ignoring the massive modes 
associated with the amplitude fluctuations, we write 

A ab (x,r, z) = u ac (x,r, z)A cd (x)vJ db (x,T, z), u e SU(2n). 

In terms of the left (i/jl) an d right (i/jr) moving fermion fields defined in the continuum limit 
around the Fermi points in the strongly coupled x-direction, the action in Eq. (f|) can be 
written as 
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5 = Ti(i) R d-i) R + i) L d + i) L ) - 2iStTr(i) R t/) L - ^l^r) 

51 = Tr(A 2 /2U) - A Tr(i> R uAu^ L + ^uAutya), 
S± = it±Ti i> R (z)i/; L (z + 1) + ii L {z)i) R {z + 1) - c.c. 

Here Tr stands for the trace over space-time as well as the replica and energy indices, 
d± = d T ±ivpd x with vp the Fermi velocity. Next, we perform a local gauge transformation, 
'4 ) lr = u vl,r, an d define the pure SU(2n) gauge fields A± = —iu^d±u = A T ± ivpA x . The 
action becomes (dropping the primes), 



S = Tr 



i/j R (d- + iA_)i) R + i) L {d + + iA+)i) L 



+ iA'Tr [ij R e iA ^/ 2+2Ae ^ L - ^Le" iA(7r/2+2Ae V 
+ Tr(A 2 /2f/) + S ± (^ LtR - uty L) .R). 



(5) 



where A' = (A 2 , + 45 2 t 2 ) 1 / 2 and 2A9 = tan- 1 (2<5t/A ). 

The final step is to integrate out the fermion fields. In order to do so, we need to 
bring the term proportional to A' in Eq. (KJ) to the standard mass term for Dirac fermions 



iA'Tr 



ipRipL — ipLipR ■ This can be done by the following chiral gauge transformation, 



^R,{L),a 



3 (-)iA aa (tt/4+ A6>) 



4>R,{L),a- 



(6) 



As a result, we encounter the well-known chiral anomaly |L2|, which arises from the Jacobian 
associated with the transformation (Bl) and leads to, 



Schirai^) = % - (| + Afl) Tre^Ad^A,, 



(7) 



where /i, v = x, r, in the transformed action. Now it is straightforward to integrate out 
the massive fermions and obtain the effective action in terms of the gauge field. Using 
the equalities Tre^Ac^A^ = (i /A)Tr e^Qd^Qd^ and Tr[v4 M ,A] 2 = Trd^Qd^Q, where Q = 
uAv) G U(2n)/U(n) x U(n), we obtain, 



S cS = ^TrdpQdpQ + ^f-Tre^QdpQduQ - ^TrQ(z + l)Q(z) - hTrAQ. (8) 



Here we have rescaled the coordinates by Ax — > x, Xvft — > r and 2^2;. The coupling 
constants cr° a have the meaning of conductivities defined on the length scale cutoff A. For the 
present network model, a° xx = (A / V^A') 2 , cr° xy = l/2 + 2A6>/vr, and a° zz = (t ± A / ^/nv F A') 2 
for small t±/vp and 5, and h = tjAo/vfUX 2 . The coupled-layers in the thermodynamic limit 
is thus described by the zero-temperature properties of the above (2+l)D quantum NLcrM. 

For t± = 0, a zz = 0. Eq. (H) reduces to N independent 2D NLcrMs discovered by 
Pruisken and coworkers for the single layer IQHE ||13|| . In this case, the term that couples to 



a xy becomes a topological quantity which produces the critical fixed points at (a xx , a xy ) = 
(const, £+1/2) for the plateau transitions, and the stable fixed points at (a xx , a xy ) = (0, £) for 
the quantum Hall states in units of e 2 /h. For the network model, the £ = 0^1 transition 
happens at 5 = where a xy = 1/2. For a xy 7^ 1/2, the conductances (a xx , a xy ) flow to (0, 1) 
for 5 > and (0, 0) for 8 < 0. At the transition, the dissipative conductance has a critical 
value a xx ~ (0.58 ± .05) p^JT5| . Notice however, for finite a° zz a xx , the system is highly 



anisotropic but three-dimensional. The cx^-term, having two derivatives, no longer contains 
nontrivial topological contributions from slowly-varying field configurations on the scale of 
the interlayer lattice spacing. 

We now present a RG study of Eq. (^). For weak interlayer tunnelings, we follow the di- 
mensional crossover analysis used in the 0(3) NLcrM description of weakly coupled quantum 



spin chains [pTql . The basic idea is that, since R = cr zz /cr xx <C 1, it is possible to consider the 
renormalization of the coupling constants in Eq. (§) in the 2D sector (x, r) independently 
until the renormalized couplings become comparable in all directions at a larger length scale 
A'. The 3D isotropic RG is switched on beyond A'. Since the scaling dimension of the Q-field 
is zero in the replica limit, this crossover takes place when Ra xx /X 2 m cr xx (\')/ A' 2 . One then 
takes the continuum limit in the ^-direction by absorbing the cutoff A' -2 into defining the 
derivatives and obtains an isotropic 3D NLaM action (plus the symmetry breaking term), 

S' cS = ^4^Trd p Qd p Q + ^^Tre^Qd.QdvQ, (9) 

o o 

with p = x,r,z and cr Q/ g(A') = a^i(X')/X r , the conductivities at cutoff A'. The important 



point is that the latter are the bare coupling constants for the subsequent 3D RG |L6 . 
Whether the system is in the insulating or metallic phase is thus determined by the renor- 
malized conductances at the end of the 2D RG. 

Since the a^-term in the continuum action Eq. (0) is no longer topological, we performed 
perturbative RG calculations to two-loop order to determine the flow of the conductance 
parameters. This approach is valid in the metallic phase where the bare conductivity is 
large for a large number of layers. For general n, we found the recursion relations for the 
conductivities under RG scale transformation, A' — > b\' 



0~ xx 



= (T xy 



l_2n I d + 2(n' + l)-—r d 

d a xx d z o± x 



(10) 
(11) 

au xx a ' u xx 

where Id = jy^' d d p/ {2ix) d [l/ (p 2 + h)], e = d — 2. Notice that the corrections to a' xy 
vanish in the replica limit n — > 0, i.e. the Hall conductivity is unrenormalized in the 3D 
regime. We will show later that this property is crucial for the quantization of the Hall 
conductance. This result should be contrasted to the one obtained in weak magnetic fields 
where the Hall conductivity is found to renormalize in the same way as the dissipative 
conductivity |T7|] . Defining the dimensionless conductances g^ v = a'X' e b e and b = e l , 
Eq. ( flO|) leads to the RG equation in the limit n — > 0, dg xx /dl = eg xx — {A/ dK d )g~^ 
Kd = 2 d ~ 1 Tc d / 2 T(d/2), consistent with the known result for the unitary NLcxM without the 



cx-cy-term in 2 + e expansions For d = 3, there is a nontrivial critical fixed point at g c = 
A/dK d e = 1/v^3tt 2 - It separates a metallic phase with g xx (\')/g c > 1 from an insulating 
phase, where g xx {X')/g c < 1- The dissipative conductivity vanishes at the metal-insulator 
transition whereas the Hall conductivity remains close to its bare value at the beginning of 
the 3D RG. The Hall conductance follows the simple Ohm's law, dg xy /dl = (d — 2)g xy . 

We now discuss the phase structure assuming the carrier density in each layer to be 
nominally the same, (i) For 5 = 0, o~° = 1/2. The individual layers are at the critical 



point for the 2D plateau transition. During the initial 2D RG, o xy does not renormalize 
and alt flows towards its finite critical value which is of order one. Thus the crossover 



length A' ~ X/y/R. For small interlayer tunneling R -C 1, A' 3> A such that cr^(A') flows 
towards the fixed point value cr^(oo) = o c xx ~ 0.55 Hl^ , |l5| . The latter is greater than the 
3D critical conductance g c derived above. Thus we conclude that when the Fermi energy 
is located at the critical point for the 2D plateau transition, an arbitrarily small interlayer 
tunneling leads to a 3D metallic state. A unique feature of the metallic phase is that the Hall 
conductivity is unrenormalized and remains close to u xy {X') down to low temperatures, (ii) 
For 5^0, the 2D RG scales towards the insulator /quantum Hall fixed points, i.e. a xx — > 0, 
whereas o 2 xy — > (8 < 0) and 1 (5 > 0) around the first plateau transition. The metallic 
phase is stable so long as a xx (X') > g c . Clearly, with decreasing (increasing) t± (\S\), at a 
critical t c ± (5 C ) where a xx (X') = g c , a metal-insulator transition takes place. For t± < t± c 
(or \5\ > 5 C ), the system is in the 3D insulating phase. To determine the phase boundary, 
notice that for 5 ^ 0, a finite localization length develops in the 2D sector, £ 2 d oc |5| _I/2d 
with u 2 d — 7/3 0,|K],|15|] . Thus the 2D RG flow stops at £ 2 d beyond which a gap would 
develop in the 2D sector. Setting the crossover length A' = one finds that the critical 
anisotropy R c ~ (A/^2d) 2 , leading to the phase boundary t c ± oc i|5|^ 2d . The width of the 
metallic phase is then given by Ws oc (t±/tY^ U2d , consistent with the numerical results of 
Chalker and Dohmen || . (iii) From the above discussion, it is clear that the metal-insulator 
transition is in the 3D unitary universality class of the Anderson transition since the Hall 
conductivity appears to be a 3D RG invariant. The two-loop RG equations imply that the 
3D localization length diverges as £3^ oc \g xx — g c \~ U3d , ^3d = l/2e. Simulations of various 
unitary models give = 1-35 ± .15 |19| which indeed agrees with the numerical value 
1.45 ± .25 obtained directly from the layered network model 0. 

We next discuss the quantization of the Hall conductance in the insulating phases. In this 
case, during the first stage of the RG, the Hall conductance a 2 xy in Eq. (|9]) flows towards the 
2D quantized values, i.e. cr X y(\') — > ie 2 /h, for large anisotropy (R 1) such that A' — > 00. 
The quantization of the 3D Hall conductance is then possible provided that a xy in Eq. (|9p 
does not renormalize in the 3D regime. Restoring the discrete sum in the z-direction, i.e., 



S 



(1/A') J dz — > X) 2 5 this term becomes S xy = Na^.y(X') Jdxdrtre^Qd^QduQ, leading to the 
Hall conductance quantization a xy —>■ iNe 2 /h. The 3D quantum Hall states are therefore 
characterized by the resistance behaviors p xx ,Pyy — > 0, p zz — > oo, and p xy = (iiV) -1 /i/e 2 . 
Due to the intervening metallic phases, the transitions between the quantized Hall plateaus 
comprise an insulator-metal and a subsequent metal-insulator transition, and have a finite 
width at low-temperatures. These results are consistent with the experimental observations 
of the IQHE in the 30-layer |§ and the more recent 200-layer GaAs/AlGaAs structures [|J. 

In contrast to 2D, where the electronic states are localized at all energies except a critical 
set of zero measure, it is the existence of metal-insulator transitions and the absence of 
localization corrections to the Hall conductivity that give rise to the quantization of the 
Hall conductance in weakly coupled layered systems. This is supported by our two-loop RG 
results in Eqs (|T0|) and (|ll|). Although these results do not form a proof to all orders in the 
perturbative RG, we believe that the evidence is sufficiently strong. 

We now briefly discuss the possible topological effects not included in the present analysis. 
The discrete (in the z-direction) action in Eq. (|§|) allows contributions from topologically- 
stable, point-singular field configurations, i. e. the hedgehogs at which the instanton number 
changes abruptly from one layer to the next |20(| . These contributions could in principle enter 
during the first stage of the RG and modify the bare parameters of the continuum action 
in Eq. @. The precise effects of the hedgehogs in the replica limit is not understood at 
the present time. Nevertheless we do not expect them to change the main results discussed 
above. 

Finally, we consider the surface states in the quantum Hall phases where the bulk lo- 
calization length is very short. The edge state supported by each layer couples together 
and forms an interesting 2D surface system decoupled from the bulk P,|2T|. In the pres- 
ence of boundaries, in addition to the bulk a xy -term, the action in Eq. (|7|) leads to an 
additional contribution (ia xy /2) § dr^trAA^, where the integral is over the boundary of the 
sample at (x = 0, L) while keeping the periodic boundary condition in r. It is easy to show 
that this surface term can be written in terms of Q(u,t,z), a smooth homotopy between 
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Sj = ^ f 1 duTiQd u Qd T Q + ?f-Trd z Qd z Q. (12) 

4 Jo o 



Q(u = 0) = Q(x = 0) and Q(u = 1) = Q(x = L). The action on the surface is then, 

'1 

4 7o 

This action is identical to the coherent state path integral action of an SU(2n) ferromagnetic 
Heisenberg spin chain with spin S = u xy /2 and exchange J = — <7 zz /<J^ y - The first term in 
Eq. (|l2j ) corresponds to the Berry phase term. The spin quantization in this case results 
from the Hall conductance quantization in the quantum Hall state. The equations of motion 
of S s i give the exactly known one-magnon dispersion valid for all n, i.e. iq T = \J\Sq%. By 
the analogy between the spin-spin correlation function and the edge electron two-particle 



Green's function [10|, it can be shown that this mode corresponds to the anisotropic diffusive 



mode on the surface, iup = —iq T + {ci Q zz /2a X y)q 2 z . Following Wegner f2~2"fl , the latter leads to 



the conductivities on the surface: af z = a° zz /2a xy and crf T (u) oc i/uj. The present approach 
to the chiral surface state is complimentary to those formulated using supersymmetric fields 
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